Frequency response in surface-potential driven electro-hydrodynamics 
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Using a Fourier approach we offer a general solution to calculations of slip velocity within the 
circuit description of the electro-hydrodynamics in a binary electrolyte confined by a plane surface 
with a modulated surface potential. We consider the case with a spatially constant intrinsic surface 
capacitance where the net flow rate is in general zero while harmonic rolls as well as time-averaged 
vortex-like components may exist depending on the spatial symmetry and extension of the surface 
potential. In general the system displays a resonance behavior at a frequency corresponding to 
the inverse RC time of the system. Different surface potentials share the common feature that 
the resonance frequency is inversely proportional to the characteristic length scale of the surface 
potential. For the asymptotic frequency dependence above resonance we find a w"^ power law for 
surface potentials with either an even or an odd symmetry. Below resonance we also find a power 
law uj" with a being positive and dependent of the properties of the surface potential. Comparing 
a tanh potential and a sech potential we qualitatively find the same slip velocity, but for the below- 
resonance frequency response the two potentials display different power law asymptotics with a = 1 
and a ~ 2, respectively. 

PACS numbers: 47.65.-|-a, 47.32.-y, 47.70.-n, 85.90.-|-li 



o 



> 

o 
m 
o 
\o 
o 



o 



I. INTRODUCTION 

The ability to manipulate liquids on niicrofluidic chips 
is essential to the functionality of micro total analysis sys- 
tems and recently there has been quite some interest 
in utilizing AC electrical surface potentials for pumping, 
manipulating, and mixing electrolytes. For an overview 
of electro-hydrodynamics and AC electro-osmosis we re- 
fer to Refs.E H II H H i, i and references therein. 

In this paper we consider the RC model in Ref. of 
an electrolyte confined by an insulating surface with an 
applied external potential of the general form 



Kxt(y,<) = Vo/(?/) exp(icjt). 



Here, lu is the driving frequency and f{y) is a dimension- 
less function representing the amplitude variations along 
the surface. For simplicity we consider a binary elec- 
trolyte with permeability e, viscosity rj, and conductivity 
a. Assuming a low Peclet number we may neglect convec- 
tion so that the electrodynamics can formally be solved 
independently of the hydrodynamics. On the other hand 
the hydrodynamics of course still depends on the electro- 
dynamics through the body force which effectively may 
be taken into account through a finite slip velocity, see 
e.g. Refs.imiaill^ and references therein. 



length scale of the Debye screening length A d , the driv- 
ing frequency is small compared to the Debye frequency 
LOjj, and the electrostatic energy is small compared to 
the thermal energy fc^T, i.e. the Debye-Hiickel approx- 
imation is valid. As shown by many, the model may be 
justified in detail starting from a non-equilibrium electro- 
hydrodynamic continuum description. For more details 
we refer to Refs. 0, ^ as well as 0, |^ and references 
therein. 

In the bulk of the liquid there is charge neutrality and 
the electrical potential 4> thus fulfills the Laplace equation 
= subject to the relation between surface charge 
ao and the potential drop across the capacitive element 

(1) 0(0, y, t) — K;xt(y, t) = croiy-, t)/Co as well as a continuity 

equation dt(TD{v,t) = -Jx{^,y,t) = cFdx4>{x = Q,y,t) re- 
lating the change in surface charge to the electrical Ohmic 



IS a series ca- 



current J. Above, Cq = (C^ 
pacitance with C's being the intrinsic surface capacitance 
and Cd — e/Az? being the Debye layer capacitance. 

We formally solve this problem by introducing the in- 
verse Fourier transform ^(Q) of /(y), see Eq. 



^(Q) 



dy' 
^/2^ 



f{y')e 



-iQy' 



(2) 



II. ELECTRODYNAMICS 

For the electrodynamics we follow Ajdari Q with the 
surface and the Debye layer each represented by a ca- 
pacitive element while the bulk-liquid is represented by 
an ohmic element. The model applies to the situation 
where spatial variations along the surface are slow on the 



From the solution of the spatially harmonic problem Q 
we may write the general solutions for the potential and 
the surface charge as 



cj){x,y,t) = Vbexp(zci;i) 

dQ -^{Q) e^p{iQy) exp(- 



■\Q\x) 



/2Tr 



1-i\Q\Ad 



(3) 
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and 

'^oiy, t) = exp{iujt) 

id 

dQ_ ,^{Q)\Q\exp{tQy) 
J-ooV2^ 1-i\Q\Ad ■ ^ ' 

Above, we have introduced the frequency dependent ef- 
fective length scale Ad given by 

AD = il + 5)^\D, (5) 

CO 

where 6 = ^ff- is the capacitance ratio and lod = ^ \s 
the Debye frequency. 



be seen by substituting the general expression for JF{Q). 
Spatial asymmetry in the applied potential is thus not 
a sufficient condition to achieve pumping and as also in- 
dicated in Ref. '3, asymmetry in the surface capacitance 
is required. Formal studies of the model for dyCg ^ 
will be an interesting subject for future studies, but here 
we will follow the lines were Cg is constant along the 
surface. The case of f{y) = cos{qy) was first studied in 
Ref.H The Fourier transform is then a simple sum of 
two Dirac delta functions, ^{Q) oc S{Q — q) + S{Q + q), 
and substituting into Eqs. {Tj) or ||SJ| we easily arrive at 
results fully consistent with those reported previously. 

IV. FREQUENCY DYNAMICS 



III. THE SLIP VELOCITY 



The slip velocity is given by Q 



(6) 



where it is implicit that real values of ctd and (j) must be 
taken. From Eqs. (O and Q we now formally get 



X Re<^ 



X Re<^ ie" 



rj oj 

'"^ dQ ^(Q)|Q|exp(^Qy) ] 
ooV2^ l-i\Q\AD J 

=° dQ' ^{Q')Q'exp{iQ'y)' 



2tt 1-z|Q'|Ad 



(7) 



and for the time- averaged slip- velocity (usiip)^ we use 
that (Re{/e^'^*}Re{5e*"*})j = |Re{/*.g} as wefi as 
^*{Q) = ^{-Q) so that 



{vs\ip{y,t)) 




dQ 



dQ' 



2?? I J-oc \/2tt 7_oo \/27r 

^(Q)|Q|exp(^Qy) ^(g')g' exp(iQ'y) 



1-i\Q\Ad 



1 + i\Q'\Ad 



(8) 



Eqs. O and © are our general results for the surface- 
potential induced hydrodynamics. Velocity and pres- 
sure fields may be obtained by solving the Navier-Stokes 
equation without body-forces subject to a slip-velocity 
condition. 

Before turning our attention to specific potentials 
we may already extract some general properties of the 
model. From Eq. © it is clear that (wsiip(2/, i))^ is a 
spatially dependent function of finite magnitude. This in 
turn also means that the external potential can induce a 
finite time-averaged velocity field in the fiuid. However, 
it is also clear that the induced net flow rate will be zero 
since {vs\ip{y,t)')^ has no spatial dc-component, as may 



In order to draw some general conclusions about 
the frequency dynamics of the hydrodynamics, apart 
from the frequency doubling jj, we imagine that the 
surface potential has a characteristic length scale 1/q. 
Since we have made an RC description of the electro- 
hydrodynamics it follows quite naturally that the system 
will have a resonance behavior at a frequency comparable 
to the inverse RC time of the system. Above resonance 
the system has an asymptotic power-law frequency re- 
sponse as we will illustrate below for a surface potential 
of even symmetry. In that case ^{Q) is a real and even 
function of Q, so that Eq. ||SJ| with Q = sq becomes 



ds ^{sq)s cos{sqy) 



V2tt %~ IS 
ds' i^^{s'q)s'sm{s'qy)' 



(9) 



where we have introduced a characteristic resonance fre- 
quency Co as well as a characteristic flow velocity u. 




(10) 



^ {s - s')n-^ + o (n-^) (11) 

we expect with fl — uj/lu that 

(wsiip)t oc w > (12) 

Likewise, for surface potentials of odd symmetry we ar- 
rive at the same conclusion. Making a similar expansion 
for the below-resonance asymptote one might speculate 
that (wsiip)j oc uj" with a power a — 2. However, with 
such an expansion the integrals are in general no longer 
convergent which in turn indicates that the power a is 
below 2 depending on the actual surface potential under 
consideration. In the following we will illustrate this by 
the aid of two explicit examples. 
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FIG. 1: Dynamics for tanh(gy) and sech(gj/) surface potentials, panels (a)-(b) and (c)-(d), respectively. Panels (a) and (c) 
show the time-averaged slip velocities at resonance while panels (b) and (d) show the corresponding frequency response. In (a) 
and (c) the solid lines show numerical evaluations of Eqs. I13I I and 1151 1. respectively, while dashed lines show the slope dy{vs Up)_^ 
at y — calculated from Eqs. 11141 and ll(ill . respectively. In (b) and (d) the solid lines show numerical evaluations of Eqs. 1141 
and II16II . respectively, while dashed lines are included to indicate asymptotics. 



A. Two-electrode configuration 

Recently, a circulating flow was observed experimen- 
tally using two adjacent electrodes and an AC field 
and in Ref. ^3 problem was studied theoretically 
by assuming that the two electrodes are infinitely close 
corresponding to a step-like model potential. However, 
such a potential is somewhat unphysical since the electric 
field diverges and instead one should take the electrode 
gap into account 14]. In the following we model such an 
electrode configuration with a smooth external potential 
of the form f{y) — taiih{qy) which incorporates the built- 
in electrode gap in terms of the length scale 1/q. Fourier 
transforming and substituting into Eq. |(SJ) we get 

/ f ,^\ -T, / n , csch(fs)ssin(sg?/) 
[Jo ^~ ~^ '^^ 

where the resonance frequency is given by tj* = tD/2. In 
order to study the maximal value of (wsiip(2/, as a 
function of frequency we note that when (wsUpCy, is 
maximal so is its slope at y = 0. If we define V{lo) = 



{qv) ^9j,(wsiip(y,<))Jy=o 

9. 9 , ("^ csch(^s)s" , ^ 

r = -2^\JI ~ J,J,), J„ = ds (14) 

where we have introduced — uj/uj* . Panel (a) of Fig.Q 
shows the time-averaged slip velocity at resonance eval- 
uated numerically from Eq. (|13|) and its slope at y = 0, 
based on a numerical evaluation of Eq. H14() . is also in- 
cluded. Panel (b) shows r(w) as a function of frequency. 
The plot confirms Eq. H12|l and suggests a linear depen- 
dence (fsiip)j c>c uj for the below-resonance frequency 
response. This may also confirmed by an asymptotic 
analysis by first noting that Ji — 1/{2Q) + 0{Q^). For 
J2 we note that dJ2/dn = -217 ds csch{^s)s^ /[fl'^ + 
4s2]2 -An/n dss/[n^ + As^ = -l/27rrj so that 
J2 ~ const. — ln(ri)/27r. For J3 we similarly note that 
dJs/dn ~ -1/8 so that J3 = const. - f]/8. When col- 
lecting terms we thus get a w-dependence to lowest order. 

B. Single-electrode configuration 

As a second example we consider single-electrode 
driven electro-hydrodynamics and as a model potential 
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we use /(j/) = sech(gy) to mimic the potential of a single 
electrode with a characteristic width 1/q. Fourier trans- 
forming and substituting into Eq. we get 



1. 



ds 



i-^ sech(|s)scos(s gy) 



' , sech(f s')s'sin(s'gy) 1 

713 h (15) 



where in this case u* = u) is the resonance frequency. In 
analogy with the previous example we get 



In = 



sech(fs)s" , ^ 



Panel (c) in Fig. ^ shows the time-averaged slip velocity 
at resonance evaluated numerically from Eq. (|15l) and the 
slope at y = based on a numerical evaluation of Eq. 1)16(1 
is also included. Panel (d) shows the frequency response 
confirming Eq. (|12|l for the above-resonance response 
while for the below-resonance response (Dsiip)^ cx; 
with a power a ~ 2. For the asymptotic analysis we 
note that = 1 -h 0{Vt^) and h = SC/tt^ + 0{9F-) 
where C ~ 0.915966 is Catalan's constant. For /i it fol- 
lows that dh/dn = -2rj/p°°dssech(fs)s/[f72 + s^Y - 
-2f7 dss/[Vi^ + s^Y = -Vl-^ so that Ii = const. - 
ln(r2) + 0{VP'). There will thus be a logarithmic correc- 
tion Lo^ ln(a>) to the uj^ dependence. 



of Ref. '3. In the case with a spatially constant intrinsic 
surface capacitance we generally prove that the net flow 
rate is zero. A finite net flow rate will require a spa- 
tially inhomogeneous capacitance so that the charge and 
electrical field dynamics are mutually out of phase. 

For the frequency dynamics we observe a resonance 
behavior at a frequency corresponding to the inverse RC 
time of the system and different surface potentials share 
the common feature that the resonance frequency is in- 
versely proportional to the characteristic length scale of 
the surface potential. Above resonance we in general find 
a uj^'^ power law for surface potentials with either even 
or odd symmetry. Below resonance we likewise find a 
power-law dependence with a being positive and de- 
pendent of the properties of the surface potential. As 
two examples we have compared tanh (odd symmetry) 
and sech (even symmetry) surface potential. Unexpect- 
edly, we qualitatively find the same slip velocity, while 
the frequency dynamics is very different for the below- 
resonance frequency response where power law asymp- 
totics with a — 1 and a ~ 2, respectively, are found. 

Our study illustrates how the frequency dynamics car- 
ries strong fingerprints of the driving surface potential 
and we believe that detailed analysis of the frequency dy- 
namics would be a way to extract more information from 
experimental data. Recent related work 0| also sug- 
gests the possibility for a strong frequency dependence 
in ac electroosmosis also in the flow topology. 



V. DISCUSSION AND CONCLUSION 

In this work we have used a Fourier transform approach 
to make a general analysis of the frequency dynamics of 
AC-induced electro-hydrodynamics within the RC model 
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